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(54) Method for determining propagation delays and method for detecting processor-clock drift 
in a communications system 



(57) A method of determining propagation-delay 
changes in the send-and-receive transmission paths 
connecting two communications units (A, B) in a com- 
munications system employs a multi-stage process in 
which, in first stage - in which it assumed that the prop- 
agation delays are equal - approximations (tp1*, tp2*) 
of the propagation delays (tp1, tp2) are derived in the 
two units, a value (dt) of the phase difference between 
the two clocks driving the respective processors of the 
two units is derived and initial values of the propagation 
delays are calculated using the approximations and the 
phase-difference value. Subsequently, in one or more 
further stages further approximations (tp 1 *, tp2*) of the 
propagation delays (tp1, tp2) are derived and, using 
these further approximations and the phase-difference 
value (dt) , further values of propagation delay are cal- 



culated, these being taken as the actual new propaga- 
tion-delay values. This process continues in a cyclical 
fashion, with changes of propagation delay being de- 
tectable by simple comparison of successive corre- 
sponding delay values. Also disclosed is a method of 
detecting drift between the two processor clocks, in 
which the method just described is basically followed 
and the propagation-delay approximations are moni- 
tored to detect opposite changes in the pair of approxi- 
mations. When opposite changes are detected, the rate 
of change is determined and, if it is below a predeter- 
mined rate, it is concluded that clock drift has occurred. 
In a preferred embodiment, once drift has been estab- 
lished an actual change measured between successive 
approximations is used to adjust the initial value of the 
clock phase difference and thereafter each new updated 
value, thereby compensating for the drift. 
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Description 

[0001] The invention relates to a method for determining the values of propagation delay in a bi-directional commu- 
nications system, and to a method for detecting drift between the independent processor clocks used in such a system. 
s The invention finds particular application in the field of power-line current and voltage measurement using relay appa- 

[0002] Power-line relay arrangements are known in which two relays - assumed designated relay A and relay B (see 
Figure 1a) - measure a value of current and/or voltage at respective ends of a power line (not shown). Each relay may 
include, as shown in Figure 1b ; a processor unit 10 driven by a clock source 11, a message-transmitting means 12 

10 and a message-receiving means 13 connected at one end to the processor 10 and at the other end to respective 
antennae 14, 15, which may in practice take the form of a common antenna. Providing measurements of the various 
electrical parameters relating to a power line 16 are current and voltage sensors 17, 18, respectively, which provide 
inputs to the processor 1 0 to be processed and finally included in a data message string to be transmitted by antenna 1 4. 
[0003] The relays operate as communications apparatus in a communications system in which the messages from 

is each relay are communicated to the other relay by wireless, normally by radio or optical fibres. In the example shown 
in Figure 1a, relay A at time tA1 transmits to relay B a message consisting of a data word, being a current vector, and 
a number of other words containing time information, including the time tA1 at which it sent off the message. Relay B 
receives this message at time tB* after a propagation delay equal to tp1 and, after a delay td, transmits its own message 
at time tB3, this message likewise consisting of a current-vector word and a corresponding number of time-related 

20 words. These words are shown to include time tB3, time tA1 which relay B learned from receipt of A's previous trans- 
mission, and time td, the time interval between tB* and tB3. Finally B's message reaches A at time tA* 
[0004] Normally, in power-line measurements it is required that the two relays compute current and other electrical 
quantities based on samples taken at the same point in time. However, it can be seen from Figure la that relay A's 
sampling points may well not coincide with relay B's (tB3 lies between tA3 and tA4), so that relay A must find a way of 

25 adjusting the value of B's vector so that it has the value it would have had if it had been sampled at time tA3 or tA4. 
To achieve this, relay A computes the total propagation delay tp1 and tp2 between the two relays as: 

tp1 +tp2 = tAMA1 - td 

30 

and, on the assumption that the two propagation delays are the same, so that: 
tp1 = tp2 = tp = (tA*- tA1 - td)/2, 

35 

calculates the time tB3* as: 



tB3* = tA* - tp. 

[0005] Relay A is then able to rotate the current vector of relay B either forwards to coincide with tA4, or backwards 
to coincide with tA3 (which presumably has been stored in a memory location to allow further processing) by an angle 
corresponding to a time difference tA4-tB3* and tB3*-tA3, respectively. This process is illustrated in Figure 2 and applies 
equally by analogy to rotation by relay B of a current vector received from relay A. 

45 [0006] While this approach works well as long as the propagation delays are equal, which in practice means using 
the same transmission route in both directions, there are occasions when it will be found necessary periodically to 
switch transmission in one direction (e.g. the direction from relay Ato relay B) to a different path, while perhaps leaving 
the other path (B to A) as it was. This will result in unequal propagation delays and consequent problematic inaccuracies 
in vector adjustment in the relays. 

so [0007] The present invention aims to address this problem of unequal propagation delays by providing a way of 
determining the values of the individual propagation delays and, consequently of detecting changes therein, and also 
tackles the consequences of drift that can arise between the clocks driving the various relays. 
[0008] According to a first aspect of the invention, there is provided a method for determining the values of respective 
propagation delays (tp1 , tp2) of send and receive paths in a communications system, said system comprising first and 

55 second spatially separated communications apparatus (A, B) each comprising a message-transmitting means, a mes- 
sage-receiving means and a processing means driven by a clock source, the clock sources having a mutual phase 
relationship (dt), said method comprising the steps of: 
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in a first stage: 

(a) establishing respective first approximations (tpl *, tp2*) of said propagation delays (tpl , tp2); 

(b) under an initial assumption that said propagation delays are equal, using said first approximations to derive 
a value (dt) representative of said phase relationship; 

(c) using said first approximations (tpl *, tp2*) and said phase-relationship value (dt) to derive first values (tpl , 
tp2) of said propagation delays, said first values being equal values; 

in each of one or more successive further stages: 

(d) establishing respective further approximations (tpl*, tp2*) of said propagation delays (tpl, tp2), and 

(e) using said further approximations (tpl*, tp2*) and said phase-relationship value (dt) to derive further values 
(tpl , tp2) of said propagation delays. 

[0009] Steps (a) and (d) may utilise the following parameters: a first transmit time (tA1 ), relative to said first apparatus' 
own clock source, at which said first apparatus (A) transmits a message to said second apparatus (B); a first receive 
time (tB*), relative to said second apparatus' own clock source, at which said second apparatus (B) receives said 
message; a second transmit time (tB5), relative to said second apparatus' own clock source, at which said second 
apparatus (B) transmits a further message to said first apparatus (A), and a second receive time (tA*), relative to said 
first apparatus' own clock source, at which said first apparatus (A) receives said further message. 
[0010] An approximation of the propagation delay (tpl) from said first apparatus (A) to said second apparatus (B) 
and an approximation of the propagation delay (tp2) from said second apparatus (B) to said first apparatus (A) may 
be derived as: 



where tpl*, tp2* are respective said approximations, tA1 is said first transmit time, tB5 is said second transmit time, 
tB* is said first receive time and tA* is said second receive time. The value representative of said phase relationship 
may be derived as: 

dt = 0.5 (tp2* - tpl*) 

where dt is said representative value. 

[0011] In said steps (c) and (e) a propagation delay (tpl) from said first apparatus (A) to said second apparatus (B) 
and a propagation delay (tp2) from said second apparatus (B) to said first apparatus (A) may be derived as: 



where tpl , tp2 are respective said propagation delays. 

[0012] Advantageously, successive further propagation-delay values may be compared with respective preceding 
propagation-delay values with a view to detecting changes in propagation delay. 

[0013] In accordance with a further aspect of the invention, there is provided a method for detecting a drift between 
first and second processor clocks of respective first and second communications apparatus (A, B) in a bi-directional 
communications system, said first and second apparatus exchanging message data including data representative of 
a time (1A1 , tB5) of transmission of said message data by each apparatus to the other, and of a time period (tdA, tdB) 
elapsed between receipt by said each apparatus of message data from the other apparatus and said transmission of 
message data by each apparatus, said transmission time (1A1, tB5) and said time period (tdA, tdB) being referenced 
to a nominal clock rate of said each apparatus, wherein each apparatus (A, B) on a cyclical basis: 
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(a) derives, from a time (tA*; tB*), referenced to its own nominal clock rate, of receipt of message data from the 
other apparatus and from said data representative of a transmission time (tA1 ; tB5) and time period (tdA; tdB), a 
first value (tp1*) approximating a delay in propagation (tp1) of said message data from said first apparatus (A) to 
the second (B) and a second value (tp2*) approximating a delay in propagation (tp2) of said message data from 

s said second apparatus (B) to the first (A); 

(b) compares said first and second approximating values (tp1*, tp2*) with corresponding values previously ob- 
tained, and if corresponding said first values (tp1 *) and corresponding said second values (tp2*) have changed in 
opposite respective directions and if a rate of change is below a predetermined rate, infers that a drift between the 
clock rates of said first and second clocks has taken place. 

[0014] Preferably, each apparatus in step (b) determines that said change in said second values (tp1*, tp2*) is mo- 
notonic before inferring that a drift has taken place. 

[0015] Advantageously, each apparatus is arranged to determine values of the propagation delays between said 
apparatus and, in a further step (c), to employ said inference that a drift has taken place to correct for errors in said 

is propagation-delay values arising as a result of said drift. 

[0016] In a first cycle at the start of a communications process between said apparatus and subsequent to said step 
(a), each apparatus may derive from said approximating values (tp1*, tp2*) a value (dt) representative of a phase 
difference between said clocks, and in said first and in subsequent cycles each apparatus may derive from said ap- 
proximating values (tp1*, tp2*) and from said phase-difference representative value (dt) said propagation delay values 

20 (tpl, tp2), each said apparatus in subsequent cycles in said step (b) detecting a predetermined amount of change in 
said approximating values (tp1*, tp2*) and in said step (c) applying that amount to modify said phase-difference rep- 
resentative value (dt). 

[0017] The invention accordingly provides a method for determining the values of respective propagation delays 
(tp1, tp2) of send and receive paths in a communications system, which method takes account of the desirability of 
25 correcting the propagation delays for drift between the clock rates, the communications system comprising first and 
second spatially separated communications apparatus (A, B) each communications apparatus comprising a message- 
transmitting means, a message receiving means and a processing means driven by a clock means, the first and second 
clock means having first and second clock rates and a mutual phase relationship (dt) therebetween, the method com- 
prising the steps of: 

30 

recording a time tAlwhen the first communications apparatus (A) sends a message to the second communications 
apparatus (B) over the send path, 

recording a time tB* when the message from the first communications apparatus (A) reaches the second commu- 
nications apparatus (B), 

35 recording a time tB5 when the second communications apparatus (B) sends a response message to the first 

communications apparatus (A) over the receive path, 

recording a time tA* when the response message reaches the first communications apparatus (A), 
calculating an approximate propagation delay tp1* in the send path as tpl * = tB* - tA1 , 
calculating an approximate propagation delay tp2* in the receive path as tp2* = tA* - tB5, 
40 calculating the phase relationship dt between the clock sources as dt = 0.5 (tp2* - tpl *) if it is an initial calculation 

of dt at a start of communications between the first and second communications apparatus, 
subsequent to the initial calculation of the phase relationship dt, 

comparing the approximate propagation delay tpl* in the send path and the approximate propagation delay 
45 tp2* in the receive path with corresponding values thereof previously obtained, 

detecting that a drift has occurred between the first and second clock rates if tpl * and tp2* change monoton- 
ically in opposite directions over a large number of calculations and if the rate of change is below a predeter- 
mined rate, 

if a clock drift is detected, deriving a clock drift value ddt from the change in tpl * and tp2* and calculating a 
so new phase relationship dt new as 

dt new = dt + ddt, and 

ss calculating the propagation delays tpl , tp2 as either 

tpl = tpl* + dt new for the sending path and tp2 = tp2* - dt new for the response path if a clock drift has been 
detected, or 
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tp1 =tp1* + dt for the sending path and tp2 = tp2*-dtforthe response path if a clock drift has not been detected. 

[0018] The first and second communications apparatus (A, B) may be relays connected to respective ends of an 
electrical power line, the relays transmitting to each other a message which includes electrical-parameter data relating 
s to said line. 

[0019] It is envisaged that occasionally communications between the relays will be lost. To take account of this, it 
may be arranged that a relative drift rate of the clocks is regularly determined and, following a break in a communications 
link between said first and second communications apparatus, said first and second propagation-delay approximating 
values are derived and from these values, and a value of drift calculated from the product of the latest said determined 
10 drift rate and a duration of said break, the respective propagation-delay values are derived. 

[0020] Embodiments of the invention will now be described, purely by way of example, with reference to the drawings, 
of which: 



Figure 1a shows a bi-directional communications system involving two power-line current relays and Figure 1b is 
is a block diagram illustrating the functional aspects of each of the relays shown in Figure 1a; 

Figure 2 is a diagram depicting the process of vector rotation carried out in each of the relays of Figure 1 ; 

Figures 3 to 6 are timing diagrams showing the process of intercommunication between two relays and illustrating 

a method according to the invention for evaluating the propagation delays therebetween, the propagation delays 

being of various values in the figures; 
20 Figure 7 is a timing diagram analogous to those shown in Figures 3 to 6 and illustrating a method according to the 

invention for detecting a negative drift occurring between the processor clocks of the two relays; 

Figure 8 is a timing diagram illustrating, in an exaggerated form, a method according to the invention for correcting 

for clock drift between the two relays; 

Figure 9 is an extended timing diagram showing the subdivision of a clock-drift correction period into four zones 
25 in order to establish a drift-correction criterion; 

Figure 10 is a diagram similar to that of Figure 8 and illustrating the clock-drift compensation method according to 
the invention in a more practical setting; and 

Figure 11 is a timing diagram similar to that in Figure 7, illustrating a method according to the invention for detecting 
a positive drift occurring between the processor clocks of the two relays. 

30 

[0021] In summary, the algorithms used in this invention are intended to enable determination of the true values of 
propagation delays in send and receive paths of communications systems as described below, even when one or more 
of the following problems are present: 

35 (1) the propagation delay times of the sending path and the receiving path of a communication channel have 

different values; 

(2) the sending path and the receiving path are of different lengths; 

(3) clock drives at each end of the communication channel are not synchronised (for example, one clock may start 
before the clock at the other end is switched on); 

40 (4) the clock drives drift apart at a certain drift rate and the pattern of drift is nonlinear; 

(5) there is a break in the communications link between the clock drives. 



[0022] A first aspect of the invention can be illustrated with reference to Figure 3. In Figure 3 it is assumed that the 
two relays A and B both sample the current at their respective ends of a power line at 5ms intervals, but that the 

45 sampling instants are not coincident - rather, relay B samples 2.5ms later (or earlier) than relay A. However, the two 
sampling rates remain identical and do not drift. Relay A sends out a message to relay B at time tA1 and this message 
is received by B at time tB*. Similarly relay B sends out its message at time tB5 and this then arrives at relay A at time 
tA*. Note that tB* and tA* are reckoned according to the individual time frames of relays B and A, respectively. The 
actual propagation delay tpl is the delay from A to B and amounts to three sampling divisions in the diagram, i.e. 15ms, 

50 and likewise delay tp2. 

[0023] However, these delays can each be broken down into two components: what might be termed a "pseudo 
propagation delay" tpl *, tp2*, which is the delay between the time of transmission of the message according to the 
transmitting relay's time frame and the time of receipt of the message by the receiving relay according to the latter's 
time frame, and a clock-difference value, dt, which represents the phase difference between the clocks of relays A and 
55 B. Thus, from the diagram, tpl * is the time between tB* and tA1, whereas tp2* is the time between tA* and tB5, dt 
being the time difference between, say, tAO* and tBO. The pseudo propagation delays are, in effect, approximations 
of the actual propagation delays. 

[0024] From Figure 3 the following equation can be derived: 
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tp2* + tp1 =tp1* + dt + dt+tp2 



dt= (tp2* + tp1 -tp1*-tp2)/2. 



[0025] It is assumed that, during the start-up period, where communications are first established, the propagation 
delays tp1 and tp2 are equal. Hence the above equation simplifies to: 



dt = (tp2*-tp1*)/2 



is [0026] To establish the value of dt it is therefore necessary to find the values of tpl * and tp2*. This is achieved by 
the following relationships: 



35 therefore, 



tpl* = tB* - tAl 

= 17.5ms - 5ms 
= 12.5ms 



tp2* = tA* - tB5 

= 42.5ms - 25ms 
= 17.5ms 



dt = (17.5 - 12.5)/2 
= 2.5ms. 



[0027] Having found dt, the actual propagation delays tpl and tp2 can be obtained as follows: 



tpl = tpl* + dt 
= 15ms 



tp2 = tp2* - dt 
= 15ms. 



55 tpl and tp2 are used by relays B and A, respectively for the purposes of time-aligning the current vector data received 
from relays A and B, respectively, as described earlier. Thus, at the very start of the communications process each 
relay calculates and records tpl * and tp2* and from these quantities dt, and then finally calculates its respective true 
propagation delay tpl , tp2. 
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[0028] Thesituation in which one transmission linkchangestoalongerpath is depicted in Figure4. Here it is assumed 
that the path delay from B to A has increased by 5ms to 20ms, while that from A to B remains unchanged. In this case, 
dt having already been calculated at the start of the communications process, that value of clock difference (dt) is 
preserved, but new values of pseudo propagation delay are calculated, and from that new values of actual propagation 
delay can be derived. 

[0029] Taking the timing details of Figure 4, then, we have: 



tpl* - tB*-tAl 

= 17.5ms - 5ms 
= 12.5ms 



tp2* = tA* - tB5 

= 47.5ms - 25ms 
= 22.5ms 

[0030] From these values, and remembering that dt = 2.5ms, we derive: 



tpl = tpl* + dt 
= 15ms 



tp2 = tp2* - dt 
= 20ms 



which are the correct values, as has been ascertained by inspection. 

[0031] A case in which each channel path increases to 20ms is shown in Figure 5. By analogy with the last example, 
we derive the pseudo propagation delays as: 



tpl* = tB* - tAl 

= 22.5ms - 5ms 
= 17.5ms 



tp2* - tA* - tB5 

= 47.5ms - 25ms 
= 22.5ms 
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tpl = tpl* + dt 
= 20ms 

tp2 = tp2* - dt 
= 20ms 

which again is the correct answer. 

[0032] In Figure 6 a situation is shown in which the two channel paths change in opposite directions: that from A to 
B increases to 20ms, while that from B to A decreases to 10ms. Applying the same equations as before: 



tpl* 



= tB*-tAl 
= 22.5s - 5ms 
= 17.5ms 



tp2* 



tA* -tB5 
37.5ms - 25ms 
12.5ms 



[0033] Hence: 

30 



tpl = tpl* + dt 

= 20ms 

35 

tp2 = tp2* - dt 

- 10ms 



which are, as before, the expected values. 

[0034] Thus, a method has been demonstrated by which the individual true propagation delays (based on the as- 
sumption that equal initial delays exist) can be calculated by the two relays and employed in deriving accurate current- 
vector rotation, even though there may have been a change in communications path in one or both directions. In 
45 practice, there will be a continuous monitoring of propagation delay on a cyclical basis, each new pair of values being 
an update on the last. By simple comparison of each update pair with the last, changes in propagation delay are readily 
detectable. 

[0035] Note that, as mentioned earlier in connection with the known relaying system, each relay relies for its calcu- 
lation of the various parameters on the receipt of the necessary time information from the other relay. Thus, as shown 
so in Figures 4, 5 and 6, the minimum information to be sent by each relay to the other, assuming that the basic clock 
difference dt has already been derived, is the value of the current vector and the time of transmission of that vector 
according to the transmitting relay's time reference. 

[0036] There remains now, however, the question as to how to solve the problem of drift in the clock rates of the 
individual relays 

55 [0037] In a practical system either or both of the clocks in the processors of the relays will tend to drift, there being 
then a net relative drift between those clocks. The situation will therefore be as shown in Figure 7, in which it is assumed 
that, purely for the sake of argument, the clock of relay A continues to sample at 5ms intervals, while that of relay B 
has decreased its sampling rate to 1 /5.5ms. This drift has been deliberately exaggerated for ease of description, a 
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more realistic drift figure being something like ±50ppm per relay. 

[0038] Although the clock of relay B is actually sampling slower than it was under the earlier examples (Figures 3-6), 
its circuitry still assumes that sampling is occurring at the same rate. Thus Figure 7 shows two time scales for relay B: 
an actual time scale, with 5.5ms sampling period, and a "time stamp" scale of 5ms period. 

[0039] As before, both relays evaluate both pseudo propagation delays over a period P, i.e. relay A, knowing tA* and 
tB5 (which is contained in B's message to A), evaluates tp2* and also, knowing tA1 and tB* (via tB5 and td B , both 
contained in B's message to A), evaluates tp1*, and similarly for relay B. 

[0040] Using the same relationships as before and comparing "time stamp" figures with "actual" figures (the latter in 
parentheses): 



tpi* = tB* - tAl 

= 15.9ms -5ms (17.5 -5ms) 
= 10.9ms (12.5ms) 



tp2* = tA* - tB5 
20 = 45ms - 25ms (45ms - 27.5ms) 

= 20ms (17.5ms) 

25 [0041] Hence, using the same value of dt as previously: 

tpl = tpl* + dt 
30 - 13.4ms (15ms) 

tp2 = tp2* - dt 
3S = 17.5ms (15ms). 



[0042] Thus it can be seen that, whereas the actual propagation delays have remained the same (15ms each), it 
appears to the relays' processing equipment as though they have changed, which could be interpreted as a change 
40 in transmission path. It is therefore necessary for the relays to be able to differentiate between a genuine change in 
propagation delay and an apparent change due to clock drift. 

[0043] In a practical embodiment of relays A and B, the internal clocks of these relays are controlled by a 16MHz 
crystal, the accuracy of which is ±50ppm. If it is assumed that the crystal of one relay is at its maximum positive error, 
while that of the other relay is at its maximum negative error, the net clock drift will be as high as 100ppm, which 
45 corresponds to a drift of 100u.s in one second. 

[0044] The minimum count resolution in the communications processor in this embodiment is 16us and therefore, 
since over a 5ms sampling period the clock will drift by a maximum of 0.5u.s, the processor will not register any change 
until 32 sampling periods have elapsed. 

[0045] Now, to detect an apparent propagation-delay change due to drift, the relay processors in the present invention 
so examine the values of tpl * and tp2* and apply a two-stage criterion test to these values, namely that if: 

(a) the values of tpl* and tp2* have changed in opposite directions, and 

(b) these changes increase steadily in the same direction over an extended period of time, 

55 then the change which has been registered will be due to clock drift and not to a change in transmission path. 

[0046] These two criteria are based on a recognition that, first of all, changes in transmission path in both directions 
and involving an increased delay in the one case and a decreased delay in the other will be comparatively rare. This 
can be seen by reference to the examples described above, where the example of "opposite changes" is only one of 
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several possibilities, one of these other possibilities involving a change in one direction only (second example) and 
another a change in both directions but in the same sense (third example), though it is appreciated that much will 
depend on the actual path selection available in a practical communications situation. 

[0047] The second criterion aims to remedy the deficiency just described in the first criterion (i.e. that there is a finite, 

s though low, probability that there will be opposite changes in path occurring more or less together) by recognising that 
a drift will be relatively slow and monotonic and cumulative in nature, whereas a path change will tend to vary in a 
random fashion and more sporadically over time. Hence the second criterion requires that a change in tp1*, tp2* in 
opposite directions occur over a large number of samples and that this change be monotonic. 
[0048] Where these criteria are satisfied, the relays then adjust the existing value of dt by a value corresponding to 

10 the change in pseudo-propagation delays at the end of the total sampling cycle, this value, however, being a multiple 
of the minimum resolution obtainable from the processors (in the example given earlier, a multiple of 16u.s), i.e. both 
relays modify their value of dt by the same amount in the same sense, positive or negative as appropriate, and dispense 
with a recalculation of tp1 and tp2. On the other hand, where these criteria are not satisfied, dt remains unmodified 
and tp1 and tp2 are recalculated to reflect a corresponding change in transmission path length. 

is [0049] This process is illustrated, purely by way of example, by reference to Figure 8. Figure 8 shows a series of 
consecutively occurring transmission routines of period P in which relay A transmits to relay B and vice-versa, there 
being a 10ms gap between receipt and subsequent transmission of a message by either relay. The actual propagation 
delays are in both cases 15ms, as in Figure 7. Also, as in Figure 7, the drift is shown in exaggerated form, i.e. a 10% 
drift (100,000ppm). The "actual" and "time-stamp" sampling points are shown on their respective time scales in the 

20 upper part of the figure together with their respective time values. It can be seen that there is a 64ms difference between 
the two after almost 15 periods P (corresponding to approximately 140 samples). 
[0050] The various relevant parameters relating to the 1 5th period P are: 

tA1 = 705ms 
25 tB* = 652ms 

tB5 = 661ms 
tA* = 745ms 

[0051] Going through the equations for the pseudo propagation delays tp1 *, tp2* as before, we have: 

30 

tpl* = tB*-tAl 
= 652 - 705 



tp2* = tA* - tB5 
•= 745-661 
= 84ms 



tpl = tpl* -I- dt 
= -53 +2.5 
= -50.5ms 
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tp2 = tp2* - dt 
= 84 - 2.5 

5 

= 81.5ms. 

[0053] The values of tp1 and tp2 are not necessarily calculated at this stage, but are shown here in order to illustrate 
10 the effect of drift on the calculations of the propagation delays and the obvious risks if such errors are not corrected. 
[0054] Now, since the relays have registered 4 counts due to the drift in the relay B clock, dt is modified by 4 x 1 6 = 
64ms in a sense such to bring the values of tp1* and tp2* closer together. 
[0055] Hence we now have: 

15 

dt' =2.5 +64 = 66. 5ms. 

[0056] If tp1 and tp2 were to be recalculated at this stage (even though it is known the delays have not changed, but 
the change in pseudo propagation delays is due to drift), we would arrive at: 

tp1 = tp1* + dt 

2S tp2 = tp2*-dt 



so tpl = -53 + 66.5 

- 13.5ms 

35 tp2 ■= 84 - 66.5 

= 17.5ms. 

[0057] The actual values are 15ms, but it can be clearly seen that the error is due to the fact that the point at which 
40 there is an exact 64ms difference between the actual and time-stamp time values in relay B is somewhere between 
period 14 and the end of period 15 in Figure 9. A similar exercise to the one just carried out but involving period 14 
would yield results for tp1 and tp2 of around 18.5ms and 12.5ms, respectively, the differences being in the opposite 
direction now. Approaching this discrepancy from a different angle, if the processor were not restricted to a minimum 
resolution of 16ms we could use a more appropriate value of time-difference to modify dt. This more appropriate value 
45 might be, for example, the difference between the time-stamp and actual values at the point corresponding to tB* in 
period 15, i.e. 717.5 - 652 = 65.5. Thus, adding now a further 1.5 to the existing 64ms fed back to change dt, we would 
arrive at: 

so 

tpl - 13.5 + 1.5 
= 15ms 

tp2 = 17.5 - 1.5 
= 16ms 
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which clearly corresponds more closely to the actual propagation delays. 

[0058] In a more preferred version of the drift-compensating method, and referring now to the more realistic drift 
rates mentioned earlier, to establish criterion (b) four "drift zones" are created in an intercommunication sequence in 
which the two relays A and B again exchange messages on a continuous cyclical basis; the zones are a 112-144 u.s 
s zone, a 144-192 u.s zone, a 192-240 u.s zone and a 240-288 u.s zone (see Figure 9). These figures in us effectively 
correspond to the increasing time differences between the "actual" and "time-stamp" time frames shown in Figures 7 
and 8, which themselves are related to the changes in pseudo propagation delay. 

[0059] To satisfy criterion (b) in this preferred arrangement, the locus of clock drift must pass through all the zones 
from the innermost zone (zone 1 ) to the outermost zone (zone 4) and must stay in each zone for a specified period, 

10 for example 40ms, before it moves on to the next zone. As soon as the "40 ms staying" condition is met for a particular 
zone, the processor allocates to the locus a flag associated with that particular zone. When the locus has collected 
four flags, one from each zone, a "clock-drift" event is registered and the before-mentioned adjustment in the value of 
dt takes place. If the drift locus passes through three zones, but fails to move on to the outermost zone and drifts back 
to the inner zones, a "clock-drift" event is not registered and the value of dt remains unmodified. If the changes of tp1* 

is or tp2* jump from the innermost to the outermost zone or beyond, the change is deemed to result from a change in 
communication path rather than from a clock drift. 

[0060] The adjustment in the value of dt in this version is illustrated with reference to Figure 10. Figure 10 corresponds 
in all essential details to Figure B, except that the drift is shown in u.s and there is a 15ms gap between receipt and 
subsequent transmission of a message by relay A and an equivalent gap for relay B of 5ms (the drawing is not shown 
20 strictly to scale). The actual propagation delays are in both cases 15ms, as in Figures 7 and 8. It can be seen that 
there is a 256us difference between the two after 50 periods P (corresponding to about 1024 current samples). 
[0061] The various relevant parameters relating to period 51 are: 

tA1 = 2548ms 
25 tB* = 2560.5ms + 256.05u.s 

tB5 = 2565.5ms + 256.55us 
tA* = 2583ms 

[0062] Going through the equations for the pseudo propagation delays tp1 *, tp2* as before, we have: 

30 

tpl* = tB* - tAl 

= 2560.756 - 2548 

35 

= 12.756 ms 



tp2*= tA*-tB5 

= 2583 - 2565.756 
= 17.244 ms 

45 Hence, 



tpl = tpl* + dt 
50 = 12.756 + 2.5 

= 15.256 ms 

55 
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tp2 = tp2* - dt 
= 17.244 - 2.5 

5 = 14.744 ms 

[0063] Again, as in Figure 8, the values of tp1 and tp2 are not necessarily calculated at this stage, but have been 
shown here in order to illustrate the effect of drift on the calculations of the propagation delays. 
w [0064] During period 51 in Figure 10, the accumulated error is 256u.s, corresponding to a 4.6 (electrical degrees - 
see Figure 2) error for a 50Hz system and a 5.53 error for a 60Hz system. 

[0065] Now, since the drift locus has passed through all four zones (periods 50 and 51 are part of the fourth zone 
since their time differences are 251 and 256 us, respectively), spending > 40ms in each zone, and a clock-drift value 
of 256 u.s has been registered (this being a multiple of the minimum processor count resolution of 16u.s), this registered 
is value is subtracted from dt to bring the values of tp1 , tp2 back to the real propagation delay times. 
[0066] Hence we now have: 

dt' = 2.5 + (-0.256) = 2.244ms. 

[0067] If tp1 and tp2 were to be recalculated at this stage (even though, again, it is known that the delays have not 
changed, but the change in pseudo propagation delays is due to drift), we would arrive at: 

2S tp1 =tpr + dt' = 12.756 + 2.244= 15ms 



tp2 = tp2*-df = 17.244-2.244= 15ms. 

30 [0068] The results of tp1 and tp2 would now be correct. 

[0069] Again, just as the earlier-described method for calculating the propagation delays was said to be cyclical, so 
also is this aspect of the invention, with drift compensation taking place on a continual basis. 
[0070] This method, then, does not correct for clock drift by changing the clock frequency of one or both relays, but 
provides compensation through software by periodically updating the value of dt arrived at initially at the very start of 

35 the communications process. Each new update acts on the value of dt as amended at the last update. 

[0071] The Table below summarises the calculation procedure performed by the algorithms described so far. 



START 


A. 


Determine basic timings: 




A1 . Record the time when End A sends a message to End B as tA1 

A2. Record the time when the End A message reaches End B as tB* 

A3. Record the time when End B sends a response message to End A as tB5 

A4. Record the time when the End B message reaches End A as tA* 


B. 


Calculate pseudo propagation delays: 




B1 . tp1 * = tB* - tA1 (sending path) 
B2. tp2* = tA* - tB5 (receiving path) 


C. 


Calculate clock drift (offset) dt: 




C1 . Ifthis is the first calculation, then dt = 0.5 (tp2* - tp1 *) 

C2. If this is not the first calculation, then dt = previous dt + clock drift 


D. 


Calculate actual propagation delays: 




D1. tp1 = tp1* + dt (sending path) 
D2. tp2 = tp2* - dt (receiving path) 


END 
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[0072] There now follows a positive example of clock drift compensation, presented in the context of a practical 
system. This is a companion example to the negative example presented with reference to Figure 7. The results from 
Figure 3 will be used for the first calculation. In particular, the clock offset dt = 2.5 ms will be used as a reference, that 
is, as "old clock offset" (dt). The maximum clock drift in this system is 40 parts per million (40 ppm), that is, the maximum 

s drift of each clock is ±20 ppm. 

[0073] The situation is shown in Figure 1 1 . Note that in Figure 7, the two time axes have different scales. In Figure 
1 1 , the two time axes have almost the same scales. This will be explained below, with the calculation process. 
[0074] In this example, Clock A is counting faster than Clock B by a rate of 40 ppm. Before End A sends out its first 
message to End B at tA1 , there is an accumulated clock drift of 2.0 ms in addition to the original clock offset dt = 2.5 

10 ms. (The 20 ms clock drift is exaggerated in order to make the diagram clearer. In the practical system, the clock drift 
is limited to no more than a few hundred microseconds, say between 250 - 500 us.) 

[0075] In Figure 3 at tA1 , the corresponding time at Clock B is 2.5 ms. In Figure 11 at tA1 , the corresponding time 
at Clock B is 4.5 ms because the drift is equivalent to a further delayed counting of Clock B. Thus: 
new clock offset, 

15 

dt = previous dt + clock drift = 2.5 ms + 2.0 ms = 4.5 ms 
[0076] The message from at End A is sent at 

tA1 = 5 ms 

and arrives at End B at 

25 

tB* = 1 5.5 ms + 0.62 u.s = 1 5.5 ms, 

where 0.62 [is is the additional clock drift between 0 ms and 15.5 ms on top of the original 2 ms. 
30 [0077] Since the drift is so small, the two time axes are almost in the same scale. The message from End B is sent at 

tB5 = 25 ms + 1.0 jlls ~ 25 ms 

35 and arrives at End A at 



tA* = 44.5 ms 

where 1 .0 u.s is the additional clock drift between 0 ms and 25 ms at End B. 
Thus, the pseudo delay 

tp1* = tB*-tA1 = 15,5-5 = 10.5 ms 



tp2* = tA* - tB5 = 44.5 - 25 = 1 9.5 ms 

[0078] It can be seen that tp1* here (10.5 ms) is smaller than tp1* derived from Figure 3 (12.5 ms). It can be seen 
so also that the tp2* here (1 9.5 ms) is bigger than the tp2* derived from Figure 3(17.5 ms). 
[0079] With these two new pseudo delays, the results are: 
Actual delay 

ss tp1 = tp1* + new dt= 10.5 + 4.5= 15 ms 



tp2 = tp2* - new dt = 1 9.5 -4.5 = 1 5 ms 
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[0080] Without the clock offset compensation, however, the results will be wrong: 
Wrong delay 

tp1 = tp1* + old dt= 10.5 + 2.5= 13 ms 



tp2 = tp2*-old dt= 19.5-2,5 = 17 ms 

w [0081] The present invention can also cater for breaks in the communications links between the two relays. Here it 
is assumed that, when the system is running as normal, the drift rate is calculated once every few seconds, e.g. every 
6 seconds. When a radio link is broken, the clocks of the two relays continue to drift, it is assumed at the same rate 
as the last-calculated rate. Thus, when the link is restored, the total net drift over the break-up period is: 

15 

Drift = Breakup period/6 sec. x Drift-rate in 6 sec 

[0082] When the link is restored, tp1*and tp2*are recalculated as described above and tpl, tp2 are derived, also as 
described above, on the basis of the recalculated tp1*, tp2*and the new drift, which is the drift dt obtaining just before 

20 breakup of the link, modified by the above drift relating to the breakup period. 

[0083] While the described embodiment of the invention has assumed certain values for the various parameters, e. 
g. a sampling period of 5 ms and a delay between receipt and transmission of a message of 15 or 5 ms, other values 
are equally possible, depending on the application. The same applies to the preferred criteria of a 40ms dwell in any 
one zone and the definition of the various zones. Also, the invention lends itself to use not just in the power-line meas- 

25 urement field but to any area in which data are to be exchanged between spatially separated equipment. Indeed, the 
invention can be employed solely for the purpose of measuring propagation delays and detecting changes therein 
even without there being any concomitant exchange of other data, such as the current-vector data of the described 
embodiment. 

30 

Claims 

1. A method for determining the values of respective propagation delays (tpl, tp2) of send and receive paths in a 
communications system, said system comprising first and second spatially separated communications apparatus 

35 (A, B) each comprising a message-transmitting means, a message-receiving means and a processing means 

driven by a clock source, the clock sources having a mutual phase relationship (dt), said method being character- 
ised by the steps of: 

in a first stage: 

(a) establishing respective first approximations (tpl*, tp2*) of said propagation delays (tpl , tp2); 

(b) under an initial assumption that said propagation delays are equal, using said first approximations to 
derive a value (dt) representative of said phase relationship; 

(c) using said first approximations (tp 1 *, tp2*) and said phase-relationship value (dt) to derive first values 
45 (tpl , tp2) of said propagation delays, said first values being equal values; 

in each of one or more successive further stages: 

(d) establishing respective further approximations (tp 1 *, tp2*) of said propagation delays (tpl, tp2), and 
so (e) using said further approximations (tpl*, tp2*) and said phase-relationship value (dt) to derive further 

values (tpl , tp2) of said propagation delays. 

2. Method according to Claim 1 , in which steps (a) and (d) utilise the following parameters: 

55 - a first transmit time (tA1), relative to said first apparatus' own clock source, at which said first apparatus (A) 

transmits a message to said second apparatus (B); 

a first receive time (tB*), relative to said second apparatus' own clock source, at which said second apparatus 
(B) receives said message; 
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a second transmit time (tB5), relative to said second apparatus' own clock source, at which said second ap- 
paratus (B) transmits a further message to said first apparatus (A), and 

a second receive time (tA*) ; relative to said first apparatus' own clock source, at which said first apparatus (A) 
receives said further message. 

Method according to Claim 2, in which an approximation of the propagation delay (tpl ) from said first apparatus 

(A) to said second apparatus (B) and an approximation of the propagation delay (tp2) from said second apparatus 

(B) to said first apparatus (A) are derived as: 

tp1* = tB*-tA1 



tp2* = tA* - tB5 

where tpl* tp2* are respective said approximations, tA1 is said first transmit time, tB5 is said second transmit 
time, tB* is said first receive time and tA* is said second receive time. 

Method according to Claim 3, in which said value representative of said phase relationship is derived as: 
dt = 0.5 (tp2*-tp1*) 

where dt is said representative value. 

Method according to Claim 4, in which in said steps (c) and (e) a propagation delay (tpl ) from said first apparatus 
(A) to said second apparatus (B) and a propagation delay (tp2) from said second apparatus (B) to said first appa- 
ratus (A) are derived as: 

tpl =tp1* + dt 



tp2 = tp2* - dt 

where tpl , tp2 are respective said propagation delays. 

Method as claimed in any one of the preceding claims, in which successive further propagation-delay values are 
compared with respective preceding propagation-delay values with a view to detecting changes in propagation 
delay. 

Method for detecting a drift between first and second processor clocks of respective first and second communica- 
tions apparatus (A, B) in a bi-directional communications system, said first and second apparatus exchanging 
message data including data representative of a time (tA1, tB5) of transmission of said message data by each 
apparatus to the other, and of a time period (tdA, tdB) elapsed between receipt by said each apparatus of message 
data from the other apparatus and said transmission of message data by each apparatus, said transmission time 
(tA1 , tB5) and said time period (tdA, tdB) being referenced to a nominal clock rate of said each apparatus, char- 
acterised in that each apparatus (A, B) on a cyclical basis: 

(a) derives, from a time (tA*; tB*), referenced to its own nominal clock rate, of receipt of message data from 
the other apparatus and from said data representative of a transmission time (tA1 ; tB5) and time period (tdA, 
tdB), a first value (tpl*) approximating a delay in propagation (tpl) of said message data from said first appa- 
ratus (A) to the second (B) and a second value (tp2*) approximating a delay in propagation (tp2) of said 
message data from said second apparatus (B) to the first (A); 

(b) compares said first and second approximating values (tpl*, tp2*) with corresponding values previously 
obtained, and if corresponding said first values (tp 1 *) and corresponding said second values (tp2*) have 
changed in opposite respective directions and if a rate of change is below a predetermined rate, infers that a 
drift between the clock rates of said first and second clocks has taken place. 
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Method according to Claim 7, in which each apparatus in step (b) determines that said change in said second 
values (tp1*, tp2*) is monotonic before inferring that a drift has taken place. 

Method according to Claim B, in which each said apparatus is arranged to determine values of the propagation 
delays between said apparatus and, in a further step (c), to employ said inference that a drift has taken place to 
correct for errors in said propagation-delay values arising as a result of said drift. 

Method according to Claim 9, in which, in a first cycle at the start of a communications process between said 
apparatus and subsequent to said step (a), each apparatus derives from said approximating values (tp1*, tp2*) a 
value (dt) representative of a phase difference between said clocks, and in said first and in subsequent cycles 
each apparatus derives from said approximating values (tp1*, tp2*) and from said phase-difference representative 
value (dt) said propagation delay values (tp1, tp2), each said apparatus in subsequent cycles in said step (b) 
detecting a predetermined amount of change in said approximating values (tp1 *, tp2*) and in said step (c) applying 
that amount to modify said phase-difference representative value (dt). 

A method for determining the values of respective propagation delays (tp1, tp2) of send and receive paths in a 
communications system, said system comprising first and second spatially separated communications apparatus 
(A, B) each communications apparatus comprising a message-transmitting means, a message receiving means 
and a processing means driven by a clock means, the first and second clock means having first and second clock 
rates and a mutual phase relationship (dt) therebetween, the method being characterised by the steps of: 

(a) recording a time tA1 when the first communications apparatus (A) sends a message to the second com- 
munications apparatus (B) over the send path, 

(b) recording a time tB* when the message from the first communications apparatus (A) reaches the second 
communications apparatus (B), 

(c) recording a time tB5 when the second communications apparatus (B) sends a response message to the 
first communications apparatus (A) over the receive path, 

(d) recording a time tA* when the response message reaches the first communications apparatus (A), 

(e) calculating an approximate propagation delay tp1* in the send path as tp1* = tB* - tA1, 

(f) calculating an approximate propagation delay tp2* in the receive path as tp2* = tA* - tB5, 

(g) calculating the phase relationship dt between the clock sources as dt = 0.5 (tp2* - tpl *) if it is an initial 
calculation of dt at a start of communications between the first and second communications apparatus, 

(h) subsequent to the initial calculation ofthe phase relationship dt, 

(I) comparing the approximate propagation delay tp1* in the send path and the approximate propagation 
delay tp2* in the receive path with corresponding values thereof previously obtained, 

(II) detecting that a drift has occurred between the first and second clock rates if tpl* and tp2* change 
monotonically in opposite directions over a large number of calculations and if the rate of change is below 
a predetermined rate, 

(III) if a clock drift is detected, deriving a clock drift value ddt from the change intp1*andtp2* and calculating 
a new phase relationship dt new as dt new = dt + ddt, and 

(i) calculating the propagation delays tpl, tp2 as either 

(I) tpl = tpl * + dt new for the sending path and tp2 = tp2* - dt new for the response path if a clock drift has 
been detected, or 

(II) tpl = tpl * + dt for the sending path and tp2 = tp2* - dt for the response path if a clock drift has not 
been detected 

Method as claimed in any one of Claims 9 to 11 , in which a relative drift rate of said clocks is regularly determined 
and, following a break in a communications link between said first and second communications apparatus, said 
first and second propagation-delay approximating values are derived and from these values, and a value of drift 
calculated from the product of the latest said determined drift rate and a duration of said break, the respective 
propagation-delay values are derived. 

Method as claimed in any one of the preceding claims, in which said first and second communications apparatus 
(A, B) are relays connected to respective ends of an electrical power line, said relays transmitting to each other a 
message which includes electrical-parameter data relating to said line. 
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Fig. 1a. 
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Fig.2. 
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